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We systematically calculate quasiparticle interference (QPI) signatures for the whole phase di- 
agram of iron-based superconductors. Impurities inherent in the sample together with ordered 
phases lead to distinct features in the QPI images that are believed to be measured in spectroscopic 
imaging-scanning tunneling microscopy (SI-STM). In the spin-density wave phase the rotational 
symmetry of the electronic structure is broken, signatures of which are also seen in the coexistence 
regime with both superconducting and magnetic order. In the superconducting regime we show how 
the different scattering behavior for magnetic and non-magnetic impurities allows to verify the s" 1 
symmetry of the order parameter. The effect of possible gap minima or nodes is discussed. 
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I. INTRODUCTION 

The relation between unconventional superconduc- 
tivity and magnetism is one of the most interest- 
ing topics in condensed-matter physics. In iron- 
based superconductors^ this question is of particu- 
lar interest because in some of these systems like 
Ba(Fei_ x Co x ) 2 As 2 and SmFeAs(Oi_ x F x ) superconduc- 
tivity (SC) and metallic antifcrromagnctism (AFM) co- 
exist homogeneously^. This coexistence is characterized 
by a competition of the two ordered states which results 
in a dramatic suppression of the magnetization below the 
SC transition temperature, T c , to the extent that reen- 
trance of the non-magnetically ordered superconducting 
phase sets in at low temperatures^. 

A number of groups have argued^— that in the ferrop- 
nictides (FP), the same electrons contribute to the su- 
perconducting condensate and to the formation of the or- 
dered magnetic moment. Furthermore, it was shown^^ 
that a coexistence of the two phases occurs only for a 
large set of parameters when the superconducting order 
is of the s H -type, i.e. the superconducting gap changes 
sign between hole and electron Fermi surface pockets 
which are located around the T and (±7r,0), (0, ±7t) 
points of the Brillouin Zone (BZ) based on 1 Fe ion per 
unit cell (so-called unfolded BZ). Therefore, a large co- 
existence region in the phase diagram of FPs indirectly 
supports the -wave symmetry of the SC order in these 
systems. 

The AFM (or SDW) phase itself shows a number of in- 
teresting properties which still have not been completely 
understood. Given the electronic structure of FPs, con- 
sisting of 2 hole and 2 electron Fermi surface pockets 
located around the T and (±7r, 0), (0, ±7r) points of the 
BZi^, it is natural to assume that AFM order emerges, 
at least partly, due to near-nesting between the disper- 
sions of holes and electrons£&ii~— . Within this scenario 
the selected magnetic order is either with momentum 
Qi = (tt,0) or Q 2 = (0,tt) [Ref. which corresponds 
to ferromagnetic order along one and antiferromagnetic 
order along the other direction. The collective spin exci- 




FIG. 1: (color online) Schematic phase diagram for 
the iron-based superconductors [Ba(Fei_ :l ;Co :!; )2As2 and 
SmFeAs(Oi_ I F I )] exhibiting a regime of coexistence of spin 
density wave (SDW) and superconducting (SC) orders^. 
Regimes (I)- (IV) refers to the various parts of the phase dia- 
gram. 



tations in the ordered state reveal pronounced anisotropy 
along x and y directions which is a consequence of the 
fact that only one of the two elliptic pockets is involved 
in the formation of the AFM stated. Similarly, the elec- 
tronic structure in the AFM ordered state is anisotropic 
and the quasiparticle interference (QPI) introduced by 
scalar nonmagnetic, as well as magnetic, impurities gives 
rise to a pronounced quasi-one dimensional periodic mod- 
ulation of the real space local density of states which, 
as we argued in Refl2l[ is in agreement with imaging- 
scanning tunneling microscopy (SI-STM) experiments in 
Ca(Fe 1 _ :c Co ;c ) 2 As 2 2 2 .. 

Overall, QPI has become a powerful experimental tool 
for elucidating the nature of the many-body states in 
novel superconductors. In the presence of impurities, 
elastic scattering mixes two quasiparticle cigenstatcs with 
momenta ki and k 2 on a contour of constant energy. The 
resulting interference at wavevector q = k 2 — ki reveals a 
modulation of the local density of states (LDOS). The in- 
terference pattern in momentum space can be visualized 
by means of the SI-STM^S. In layered cuprates the analy- 
sis of the QPI has provided details of the band structure, 
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the nature of the superconducting gap, or other compet- 
ing orders 2 ^—. 

In the superconducting state of iron-based supercon- 
ductors a recent SI-STM experiment claimed to unam- 
biguously identify the order parameter symmetry to be 
of the s" 1 character—. Theoretical analyzes performed 
previously have also shown a pronounced difference in 
interference patterns characteristic of a simple s-wave 
symmetry and an extended -wave symmetry [see 
Refs03(f3c|. In particular, in the presence of vortices 
an additional source of scattering either suppresses or en- 
hances intra- and inter-pocket scattering depending on 
the symmetry of the superconducting gap2H, a feature 
found also in the experimental data3£. At the same time, 
the experimental results have been put into question and 
argued instead to arise from the Bragg scattering and not 
due to QPI22. 

In what follows we address signatures of the different 
orders on the electronic structure as seen in the distinc- 
tive features probed by QPI . In particular, we extend 
our previous analysis of the QPI in the AFM stated to 
the entire phase diagram of the FPs assuming that the 
AFM and SC arise from the same quasiparticles. We 
present the results for QPI in the coexistence regime of 
AFM+SC with s- and extended s-wave (s -1 ) symmetry 
as well as for the pure superconducting state. We fur- 
ther investigate the role of higher harmonics (non-trivial 
angular dependence along the Fermi surface) in the ex- 
tended s— wave case and, in particular, we address the 
question whether a nodeless or nodal -wave symme- 
try can be identified within SI-STM. The aim of our anal- 
ysis is to find subtle features of the various many-body 
states in iron-based superconductors that can be probed 
by STM. 

The paper is organized as follows. The theoretical cal- 
culations to obtain the local density of state (LDOS) are 
based on a four band model, and they are explained in 
Sect. II. Based on our model we numerically analyze the 
QPI signatures in the various phases in Sect. III. We 
conclude the paper by a summary (Sect. IV). 



II. THEORETICAL MODEL 

We employ an effective mean-field four band model 
with two circular hole pockets centered around the (0, 0) 
point of the unfolded BZ (a-bands) and two elliptic elec- 
tron pockets centered at (±7r, 0) and (0, ±7r) points of the 
BZ, respectively (/3-bands) 2 ^: 
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tp ([1 + e] cos(k x + 7r) + [1 - e] cos(fc y )) - 
tp ([1 — e] cos(k x ) + [1 + e] cos(fc y + n)) — H\. e accounts 
for the ellipticity of the electron pockets. Following our 
previous analysis of the spin excitations and QPI in 
the magnetically ordered state, we use Fermi velocities 
and size of the Fermi pockets based on RefsJ^, namely 
t a = OMeV, t fj = -0.68eV, fx ai = l.UeV, e = 35meV, 
/i„ 2 = lAAeV, fii = -1.23e^, and e = 0.5. For these 
values, the Fermi velocities are 0.5eVa for the ai-band, 
where a is the lattice spacing, and v x = 0.27eVa and 
v y = 0A9eVa along x- and y-directions for the /3i-band, 
and vice versa for /?2 • We use a x = a y = a = 1. 

Here, is the superconducting gap for band i. 
Symmetry and structure of the superconducting gap 
in Fe— based superconductors have been subject of nu- 
merous experimental and theoretical papers in recent 
year a 1 ' 11 ^ 2 ' 36 " — . There is a growing consensus that the 
gap has an extended s— wave symmetry - it belongs to 
a symmetric A\ g representation of the D^h symmetry 
group of a square lattice. Gap values along electron and 
hole Fermi surfaces (FS) are of opposite signs. 

A more subtle issue is whether the gap has nodes. This 
is not a symmetry issue as adding higher harmonics to 
the extended s— wave gap yields stronger momentum de- 
pendence of the latter on the FS. Quite generally one can 
write the superconducting gap in the form 

Ak = A [cosk x cos fey] + A'[cos(k x ) + cos(/c y )]/2. 

For A' = the superconducting gap can be roughly ap- 
proximated by a constant on the hole and electron FS 
pockets but with opposite signs, A e = —A' 1 ~ —A. 
Increasing A' does not significantly change the momen- 
tum dependence of the gap on the FS of the hole pock- 
ets which can be approximated as A Qi = A h while 
it introduces substantial angular dependence of the su- 
perconducting gap on the electron pockets in the form 
A^^i/j) = A e (l ± bcos2cj)). Here, <\> is the angle mea- 
sured from the line connecting the two electron FSs, i.e. 
the line between (tt, 0) and (0, it) points of the first BZ 
and b ~ 2A'/A. Such A ft (<£) has no nodes if b < 1 but 
"accidental" nodes do appear when b > 1. Their posi- 
tions are determined by cos2</> = 1/b and the latter shift 
continuously with 6, so that the node's location is not 
selected by any symmetry 

FPs are multi-orbital systems and the orbital nature of 
the quasiparticle states introduces a sizable variation of 
the interaction along the Fermi surface. This manifests 
itself as a cos 2<f> component of the interaction. In partic- 
ular, b grows upon inclusion of the intra-band Coulomb 
repulsion into the gap equation because this term couples 
to the gap averaged over the FS and hence reduces angle- 
independent gap components but does not affect cos 2<f> 
components (Ref f38h . As a consequence, b becomes pro- 
gressively larger as the system moves further away from 
the SDW phase and the effect of intra-band repulsion 
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grows. In other words ovcrdopod FPs are more likely to 
have nodes in the gap. 

We introduce the experimentally observed (it, 0) SDW 
order parameter, within a standard mean-field approx- 
imation: W oc Z) P ( C Q lP ,5 c /3ip+Qi7 < ?<57}- In this state, 
one of the a fermions couples with only one band of f3 
fermions, leaving the other hole and electron bands - 
and hence their electron and hole FSs - unaffected by 
the SDW. Without loss of generality we direct W along 
the z- quantization axis. 

The actual QPI which is believed to be measured in 
SI-STM22 arises from quasiparticle scattering by pertur- 
bations internal to the sample such as non-magnetic or 
magnetic impurities. We perform a standard analysis 
of such processes based on a T-matrix description^. In 
particular, we introduce an impurity term in the Hamil- 
tonian 
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where Vj^., and J" a , represent the non-magnetic 



and the magnetic point-like scattering between the 
electrons in the bands i, and %' respectively. In 
the following we set the quantization axis of the 
magnetic impurity along the z-direction. Following 
Ref. [33l . we define the new Nambu spinor as ^ = 

( C l 2 kt ! C aikf C /3ikf C /3 2 kt' c £*2-k|j C ai-H, C^-k^, C^ 2 -kJ.) 

where now wc measure the momentum transfer relative 
to the interpocket momentum k + Qi, e.g. c^k+Qio-- 
Therefore, the Hamiltonian can be written as 
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where by defining V^ k , = -fu ; V^., = -fUQ and J" Z S Z = 

7'uo; Jfz &z = i'uq, the matrices /3k and £/kk' are defined 
as 
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The spin index is a = ±1 for spin up and down. The 
overall a in front of the energy matrix appears due to 
the Nambu structure^. Here, we assume that the intra- 
band impurity scattering, uq differs from the intcrband 
scattering between the bands separated by a large Q, 
uq. The Green function (GF) matrix is obtained via 

GW(t) = -<TV> k (T)$ k ,(0)>, whence 

GWK) = G k (u; n )[<W +*kk'(w«)GtK)], (4) 

where G^(w„) = (iuo n — is the bare GF of the 

conduction electrons. Solving the Dyson equation for the 
T-matrix 

ikk'(w„) = ?7kk' + ^2 ^kk"Gk"(w„)tk"k'( w ™), ( 5 ) 

k" 

the LDOS is obtained via analytic continuation iuj n — > 
U) + ia0 + according to 

N c (u,r) = -hraTr[G(r,r,u n )] iuiMi(r0+ . (6) 

We recall that the interference between incident and scat- 
tered waves gives rise to the spatial modulation of the 
amplitude of the total wave. This is then reflected in the 
LDOS. We further set u = 0.3t a - 



III. RESULTS 

In what follows we discuss systematically the differ- 
ent ordered phases of the iron pnictides as shown in 
the phase diagram of Fig.([T]) and present our calcu- 
lated results for the total spectral function [Eqs.Q- 
©]. The spectral function Im Tr Go a (k, a;) of 
the clean system is always shown on the left panel. 
This allows to trace easily the most important scat- 
tering vectors that appear in the QPI corrections 
£ k(T Im Tr [G 0CT (k, u)t„ (k, k + q, u)G 0(T (k + q, u)] . We 
discuss the QPI maps for the case of non-magnetic and 
magnetic impurities and show the difference between 
them. The latter is particular important for identifying 
the symmetry of the superconducting gap. 



A. Spin-Density wave phase 

First, we review our results for the QPI in the SDW 
phase and its changes due to doping [regime (I) in Fig.([T|)] 
by showing in Figs.© and Figs.® the QPI for zero and 
7% electron doping, respectively. A consequence of SDW 
order in parent iron-based superconductors is that the 
C4 symmetry of the normal state is continuously bro- 
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FIG. 2: (color online) Constant energy intensity maps at u = 
20meV (first row) and u; = —20meV (second row) of the spec- 
tral density, ImTr Go CT (k, cj), (first panel from the left) 
and QPI, X; kCT ImTr[G 0CT (k,o;)r; (T (k,k + q,a;)G, )CT (k + q,a;)] 
for non-magnetic (second panel) and magnetic (third panel) 
impurities obtained as described in the text. The difference 
between QPI for non-magnetic and magnetic impurities is dis- 
played in the right panel. The color bars refer to the intensity 
in units of states/eV. The SDW order parameter is fixed at 
W = AQmeV for zero doping. 
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FIG. 3: (color online) Constant energy intensity maps at uj = 
20meV (first row) and w = — 20meV (second row) of the 
spectral density and QPI, for non-magnetic (second panel) 
and magnetic (third panel) impurities obtained as described 
in the text. The difference between QPI for non-magnetic and 
magnetic impurities is displayed in the right panel. For SDW 
(region I in Fig[T} with W = 40meV and with 7% e-doping, 
for uj = 20meV (first row) and w = — 20meV (second row). 
The color bars refer to the intensity in units of states/eV. The 
SDW order parameter is fixed atW = 40meV for 7% electron 
doping. 



ken. The resulting Fermi surface consists of small pockets 
which originate due to mixing of one of the hole pock- 
ets centered around the T— point and the elliptic elec- 
tron pocket centered around the (0, tt) point of the BZ 
as shown in the left panel. The scattering between the 
resulting boomerang-like structures (the most character- 



istic wave vectors qi and q2 are indicated by arrows) 
dominates the QPI at low energies and occurs for both 
non-magnetic and magnetic impurity scatterin g 21 ! 22 . 

Note that this reduction from C4 to C2 symmetry in 
the magnetically ordered state was originally interpreted 
as a sign of the underlying quasi-one-dimensional elec- 
tronic structure and the electronic nematic orde r 22 ' 42 . 
In our case, however, it appears just as a result of the 
translational symmetry breaking induced by the (0, 71") 
antifcrromagnetic state on the two-dimensional electronic 
structure. We also find that the C4 symmetry is restored 
for bias energies exceeding twice the SDW gap value. 
Recently this prediction was confirmed experimentally 4 ^ 
again signifying the dominant role of the magnetic order 
in breaking the C4 symmetry. 

For large enough W the FSs of the bands involved in 
the SDW are completely gapped and the same is true for 
the spectral densities at low energies. In this case, QPI 
will be determined by the bands which are not involved 
in the SDW, and, therefore, its structure will not show 
strong quasi-lD character. This may explain why the 
pure AF SDW state does not show any ^-symmetric 
structure in the parent compounds where the magnetic 
moment (and the corresponding SDW gap) is quite large. 
Only when it is reduced upon doping, the bands involved 
in the SDW are located close to the FS so that the QPI 
structure, described above, becomes visible. Another ap- 
parent effect is the particle-hole asymmetry of the QPI 
which is the result of a rotation of the SDW induced 
pockets from positive to negative energies by 90° degrees. 
The asymmetry is 90° rotated in the case of electron dop- 
ing (see Fig J3]) . Due to the larger size of the electron 
pockets at (±tt, 0) the SDW induced small pockets are 
pronounced and dominate the scattering features in the 
QPI maps. For low negative energies the electron pocket 
is still present and the 90° rotation of the one dimensional 
features shifts to lower energies. 



B. Superconducting phase 

Next, we discuss the QPI maps in the superconducting 
state. As mentioned above, the intensity of the QPI is af- 
fected by the relative sign of the order parameter between 
the FSs involved in the scattering via the Bogolyubov 
coherence factors representing the formation of Cooper 
pairs and new quasiparticles, which are coherent super- 
positions of electrons and holes. Coherence factors char- 
acterize how the scattering of a superconducting quasi- 
particle off a given scatterer differs from the scattering 
of a bare electron off the same scatterer. Generally, the 
momentum-dependent order parameter enters into the 
expression for the coherence factor so that studies of scat- 
terings of quasiparticles with different momenta can pro- 
vide important information on the phase of the supercon- 
ducting order parameter in momentum space. Originally, 
this idea was put forward for the cuprates^— . In par- 
ticular, for potential (non-magnetic) scattering at wave 
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FIG. 4: (color online) The spectral density and the QPI maps 
for non-magnetic and magnetic impurity scattering (from left 
to right) at 20meV in SC state with nodeless s + ~-wave (up- 
per panel) and isotropic s ++ -wave (lower panel) symmetries of 
the superconducting gap. We set A = 20meV, 13% e-doping. 
The difference between QPI for non-magnetic and magnetic 
impurities is displayed on the right panel. The arrows show 
scattering between electron and hole bands at qo = (7T, 0), 
(0, 7r). The inset in the upper panel shows the the experimen- 
tal data taken from Refl28l. 



vector q the corresponding coherence factor is smaller for 
those q which connect parts of the FS with the same sign 
of the SC gap. 

Therefore, the QPI intensity appears only for sign- 
reversing momenta. In iron-based superconductors with 
-wave symmetry of the superconducting order pa- 
rameter, this occurs for (7r, 0) and (0, 7r) momenta. Then 
for magnetic (time-reversal) scattering, the QPI is neg- 
ative for the wave vectors (ir, 0) and (0, tt) and positive 
for momenta (±7r, ±7r) as the latter connects FSs with 
the same signs of the order parameter. Overall, the ef- 
fect of the sign change of the superconducting gap can 
be most efficiently seen if one plots the difference of QPI 
for non-magnetic and magnetic impurities. 

As a result, a comparison of the QPI maps for mag- 
netic and non-magnetic scattering yields important in- 
formation on the symmetry of the superconducting order 
parameter. In practice, the most straightforward way to 
perform this comparison in a type II superconductor is 
to apply an external magnetic field where the resulting 
vortices act in part as magnetic scatterer a 28 i 45 . Thus, by 
comparing the results with and without magnetic field 
the symmetry of the superconducting order can be de- 
duced. 

In FigH]we show the QPI maps for the nodeless s H - 
wave (A' =0) and isotropic s ++ -wave symmetry. Clearly 
by subtracting the QPI for the non-magnetic and mag- 
netic scattering shown in the right panel, the difference 
between the two gaps becomes apparent with regard to 
the interband scattering (i) at q = (±7r, 0) or (0, ±7r) 
which occurs for the scattering between electron and hole 
bands and (ii) diagonal q^-scattering between two elec- 
tron pockets. As the sign of the superconducting gap is 



opposite between the a and the /3-bands the difference 
in the QPI at qo between non-magnetic and magneti- 
cally induced scattering is positive. By contrast, that 
same difference in the scattering at qi is negative as the 
order parameter on the electron pockets possesses the 
same sign. Both features arc consistent with the QPI 
measured in the superconducting state of Fe(Se,Te) and 
Ba(Fei_ a; Co a: )2As2 compound a 28 ' 47 . In particular, in the 
inset of FigHwe show the SI-STM data taken at 10T in 
FeTeo.4Seo.6 compound. Observe the difference in the 
intensities in the external magnetic field of 10T for the 
scattering at qo and q'j which is consistent with our cal- 
culations. Note also that in the difference plot the details 
of the Fermi surface scattering are very weakly visible, in- 
cluding the scattering between the electron pockets. In- 
stead the largest effect occurs due the sign change of the 
gap and the corresponding coherence factors. For com- 
parison, in FigH] (lower panel) we show the same cal- 
culations for s ++ -wave order parameter. The striking 
distinction with the s H -wave case is that the difference 
between non-magnetic and magnetic impurity induced 
QPI is, as expected, featureless. 




FIG. 5: (color online) Calculated spectral density in the 
superconducting state with nodal s +_ -wave symmetry for 
A' = 3A. The arrows refer to some intra and inter elec- 
tron pocket scattering wavevectors (i.e., wavevectors joining 
regions with large DOS). Observe that qi, and q^, q^, qi 
are the scattering wave vectors for the same and the opposite 
signs of the superconducting order parameter, respectively. 
The evolution of the banana-shape structures is shown for +5 
(black curve), +10 (brown curve), +15 (green curve), +20 
(red curve), and +30 (blue curve) meV. 

In the next step we show the effect of higher harmonics 
in the s" 1 -wave function. By taking a non-zero A', gap 
minima or even nodes form on the electron pockets. In 
this regard the evolution of the spectral density with en- 
ergy will resemble somewhat the structure found in nodal 
c?-wave superconductors. In particular, due to the nodal 
structure of the superconducting gap, the spectral den- 
sity with increasing energy will show a banana-shaped 
structure with new scattering momenta associated with 
the large density of states at the edges of the banana as 
shown in Fig [5] We selectively show only some of the 
scattering wave vectors on this plot which we later iden- 
tify on the QPI maps. Note that q^ refers to the scat- 
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FIG. 6: (color online) (color online) The spectral density 
and the QPI maps for non-magnetic and magnetic impurity 
scattering (from left to right) at +20meV (upper panel) and 
+ 10meV (lower panel) in SC state with nodal s + ~-wave order 
parameter of the superconducting gap (A' = 3A) and 13% of 
the electron doping. The difference between QPI for non- 
magnetic and magnetic impurities is displayed on the right 
panel. 



Note, to prove unambiguously whether these originate 
from the nodal s" 1 -wave order parameter measurements 
in an external magnetic field are required. 



C. Coexistence phase with SDW and s + SC order 
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tcring with the same sign of the order parameter while 
other wavevectors, q^, qjj, and q^ connect regions with 
the superconducting order parameter of opposite sign. In 
principle these additional features should be most pro- 
nounced in the difference of QPI intensities for magnetic 
and non-magnetic scattering. 

The results for the spectral density and QPI maps 
are shown in Fig|6j Due to relatively large quasiparticle 
lifetime used in the calculations (~ 4meV) the banana- 
shaped structures are not well resolved and overall the re- 
sulting QPI does not have well-defined momenta as those 
shown in Fig{6j However, the anisotropy of the supercon- 
ducting gap is clearly observed in the difference between 
magnetic and non-magnetic impurity induced QPI, as 
shown on the right panel. First, we again find the scat- 
tering along the diagonal momentum at q' 2 which was 
also present in the case of a nodeless s" 1 -wave symme- 
try. It is, however, additionally enhanced due to banana- 
shape structure and the corresponding high density of 
states at the edges of these bananas. Similarly, the inter- 
band induced scattering between the a— and the /3-bands 
(indicated by qo) acquires more structure. In addition, 
the maximum of intensity shifts from (±tt, 0) and (0, ±7t) 
to the smaller momenta due to the higher harmonics in 
the s" 1 -wave gap. Additional peaks at small momenta 
include the scattering at q^ within electron pockets and 
also the scattering between electron and hole pockets con- 
necting states with opposite sign of the superconducting 
gap shown by q5. Overall, the structure of QPI becomes 
considerably richer, which should allow identification of 
the nodal structure of s"' -wave superconducting order. 
One has to bear in mind, however, that the higher har- 
monics in the s" 1 -wave symmetry with corresponding 
scattering at q^ and q5 should appear at larger doping. 



FIG. 7: (color online) Calculated spectral density and the 
QPI maps at ui = +20meV in the coexistence regime of SDW 
+ SC. The upper panel refers to the regime II of the phase 
diagram (Fig[fj) with W = 40meV, A = 20meV and 10% 
e-doping. The lower panel describes the regime III with W = 
WmeV, A = 20meV and 12% e-doping. 

Finally, we consider the coexistence regime with both 
microscopic SDW and s' 1 -wave superconducting order. 
In Fig[7] we examine two different situations which cor- 
respond to regimes (II) and (III) of the phase diagram 
shown in FigHJ respectively. In regime (II) the SDW 
order parameter is larger than the SC one and we set 
W = 2 A. In regime (III), which corresponds to higher 
doping, the situation reverses and we use W = A/2. 
For the sake of simplicity and also because it is expected 
that the nodal -wave order occurs only at higher 
dopin g 40 ' 41 we neglect the effect of the higher harmonics 
in the gap, i.e. we put A' = 0. 

For the case of a larger SDW gap, the scattering has 
a pronounced C2 symmetry and the only effect of su- 
perconductivity shown in the QPI, as compared to the 
pure SDW case, is to reduce the overall intensity of the 
scattering peaks (note the reduced intensity on the color 
bars). The latter occurs due to additional gapping of the 
bands. The structures visible in the upper panel of FigJT] 
arc similar to those seen Figj^l with some extra smearing 
due to the superconducting gap. There are, however, ad- 
ditional structures which become particularly visible in 
the difference between QPI maps for magnetic and non- 
magnetic scattering, which are signatures of the nodeless 
-wave gap. Note that these features strengthen in 
the case when the SC gap is larger than the SDW one, 
as is seen in the lower panel of FigJTJ 

The origin of these additional structures can be traced 
back to the effect of the SC s" 1 order which is super- 



7 



imposed on the SDW order. First, notice that in the 
difference plot there is always an enhanced intensity at 
qo = (±7r, 0) momenta. This is due to the fact that 
the electron and the hole band not involved in the SDW 
state with (0, tt) order have a different sign of the su- 
perconducting gap. Then for the very same reason as 
in the pure SC state with nodeless s" 1 -wave order pa- 
rameter, the difference between the magnetic and non- 
magnetic impurity-induced QPI shows the strongest fea- 
ture for intcrband scattering at (±7r, 0) momenta. At the 
same time, new features arise due to the fact that one of 
the electron and one of the hole bands are mixed in the 
SDW state with (0, n) ordering. However, both of these 
bands still possess an opposite sign of the superconduct- 
ing gap, thus leading to a structure in the difference plot 
at small momenta. Remember that in the pure SC state 
with s" 1 -wave order, these peaks would occur for (0, ±7r) 
momentum but now this is a reciprocal wave vector of the 
SDW state. As a result, this scattering is 'translated' to 
that around q ~ momentum. This clearly shows that 
the regime of microscopic coexistence of SDW+SC can 
be nicely observed in the FS-STM data. The same occurs 
for the scattering between two electron (/3)-bands. Orig- 
inally along the diagonal of the BZ (q^ in the nodeless 
-wave symmetry), it is now shifted to a different mo- 
mentum due to the folding of the electron band involved 
in the SDW formation. 



IV. CONCLUSION 

In this paper we have systematically calculated quasi- 
particlc interference effects due to magnetic and nonmag- 
netic impurities for the whole phase diagram of the iron 
pnictide superconductors. We have shown that in the 
SDW phase the QPI shows quasi one-dimensional fea- 



tures as measured in recent experiments^. The C4 sym- 
metry is restored for energies larger than twice the SDW 
gap value. Furthermore, analyzing the superconducting 
state we have shown that the nodeless as well as the 
nodal -wave symmetry of the superconducting gap 
can be clearly identified in SI-STM experiments and dis- 
tinguished from the isotropic s-wave gap. In particular, 
the scattering between the electron and the hole pockets 
at q = (±7r, 0) [(0, ±7r)] as well as the scattering between 
the electron pockets at q^ becomes quite pronounced in 
the difference of the QPI maps between magnetic and 
non-magnetic impurities, a feature found also in the ex- 
perimental data3£. We have demonstrated that a non- 
trivial angular dependence of the s H -wave gap induced 
by higher harmonics results in banana-shape structures 
in the QPI maps not present in the simplest version of 
-wave symmetry. A large density of states associ- 
ated with the edges of these bananas is responsible for 
an additional low q-scattcring along the bond direction 
as well diagonal q-scattering patterns in the QPI. Again 
these scatterings are especially pronounced in the differ- 
ence plots between the magnetic and non-magnetic im- 
purity scattering. In the regime of coexistence of SDW 
and SC with s H -wave order our study has revealed addi- 
tional characteristic features of the QPI which may help 
to identify this phase in iron-based superconductors. It 
includes: (a) the C2 anisotropy of the QPI maps in the 
coexistence regime and (b) the low q-interband scatter- 
ing for the bands with opposite sign of the s" 1 super- 
conducting gap. 
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